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PREFACE 


This final report is the summary of a two-year effort portraying 
the development, calibration and verification of a one-dimensional variable 
cross-sectional area numerical model. This is one of the publications 
of the thermal pollution group at the University of Miami. 

This effort covers the applications of this model to Cayuga Lake 
In New York and Lake Keowee In South Carolina. The conclusions and 
derivations are accompanied by an abundance of figures and tables. 

These figures and tables are presented in such a way as to make it 
possible for use in calibrating or verifying other one-dimensional models. 

This work was conducted under funding from National Aeronautics 
and Space Administration (NASA), Kenedy Space Center and the Environ- 
mental Protection Agency (EPA). 
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ABSTRACT 


A one-dimensional model for studying the thermal dynamics of cool- 
ing lakes has been developed and verified. The model is essentially a 
set of partial differential equations which are solved by finite difference 
methods. The model includes the effects of variation of cross-sectional 
area with depth, surface heating due to solar radiation absorbed at the 
upper layer and internal heating due to the transmission of solar radia- 
tion to the sub-surface layers. The exchange of mechanical energy be- 
tween the lake and the atmosphere is included through the coupling <»f 
thermal diffusivity and wind speed. The effects of discharge and intake 
by power plants are also included. 

The numerical model was calibrated by applying it to Cayuga Lake. 
The performance was very good. The model was then verified through 
a long term simulation using Lake Keowee data base. The comparison 
between measured and predicted vertical temperature profiles for the 
nine years is good. The physical limnology of Lake Keowee is presented 
through a set of graphical representations of the measured data base. 
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SECTION 1 


INTRODUCTION 


Deep bodies of water provide a convenient source of condenser 
cooling water supply to electric generating power plants. One of the 
problems associated with this, that has gained increasing importance in 
recent years is the thermal pollution caused by the discharge of waste heat 
from these power plants. The degradation of the quality of these waters 
usually occurs through the direct influence of the Increased temperature 
on aquatic life or through the lowering of the amount of dissolved oxygen 
and other biochemical effects. 

In temperate regions most deep bodies of water develop a thermo- 
cline during their annual heating cycle. A warmer epilimnion at the top 
Is isolated from a cooler hypollmnion below by severe stable gradients 
(thermocline) . Convective transport and heat addition caused by power 
plant discharge result in disturbances in the thermocline. 

The formation time, phasing, depth and severity of the thermocline 
are crucial factors affecting the bio-chemical processes in an aquatic 
ecosystem. The nutrient levels, species spectra and physical character- 
istics are quite different in the two distinct domains below and above the 
thermocline. This stratification is believed to be caused by nonlinear 
interaction between the wind- generated turbulence and stable buoyancy 
gradients. While being heated from above, a basin forms stable strati- 
fication thereby inhibiting wind-generated turbulence. During early 
spring, most temperate lakes exhibit a nearly homothermal temperature 
distribution with a temperature of about 4*^C (which is the temperature 
of maximum density for water) extending all the way to the bottom. As 
the air above the lake begins to warm, the lake receives heat, at an 
increasingly rapid rate. During the early part of the warming season, 
the lake continues to remain nearly homothermal, since the heat that is 
received at the surface layers is transported to the deeper layers by 
wind-induced currents and turbulence. As the rate of heating of the 
lake continues to increase, the rate at which heat is received at the 
surface layers soon exceeds the rate of heat removal to the deeper lay- 
ers, and the temperature of the surface layers begins to increase. 

During this period the temperature decreases monotonically with increas- 
ing depth. As heating continues, a point of inflexion develops in the 
temperature profile separating an upper well-mixed layer from a lower 
well-mixed layer. The region around the point of inflexion is a region 
of intense stable temperature gradient and consequently, low turbulence 
levels. Turbulent diffusion through this region is minimal. Further 
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heating merely increases the temperature of the well-mixed upper layer. 

As cooling begins, the wind mixing is augmented by -convective 
mixing caused by static instability and the thermocline (the region con- 
taining the point of inflexion) recedes into the deeper layers of the lake. 
The epilimnion cools until the lake reaches its minimum heat content and 
near homothermal conditions result. 

This report describes the application of a one-dimensional model in 
predicting the above phenomena. This model was calibrated using data 
from Cayuga Lake and the calibrated model was applied to Lake Keowee. 
The purpose of the effort is to develop a predictive ability for long-term 
thermal behavior of cooling lakes. 
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SECTION 2 
CONCLUSIONS 


A one-dimensional model which includes area change with depth, 
verticai convection, varying diffusivity, thermal discharges, and internal 
absorption of radiation has been developed. The model was calibrated 
using Cayuga Lake data base and verified with Lake Keowee data base. 

The application to Cayuga Lake indicates excellent performance. 

A comparison of cylindrical and paraboloid approximations for the lake 
indicate significant differences In thermocline depth, eddy-diffusivity, 
and temperatures at mid-depths. This shows that effects of area change 
with depth are not negligible. The insignificant surface temperature 
difference is attributable to the fact that equal surface areas were used 
in both simulations. However, these effects will be more pronounced 
in real basins where decrease in area with depth is more severe than 
the linear variation for the paraboloid case. 

The long term (nine-year) simulations of temperature profiles, 
formation and decay of the thermocline for Lake Keowee compare fairly 
well with measured data. For these simulations, the model did not need 
recalibration. While the influence of discharge and ambient conditions 
were satisfactorily modeled, the calculated and measured profiles had 
much larger differences for Lake Keowee than for the one-year simulation 
in Cayuga Lake. Two reasons are believed to cause this relatively great- 
er inaccuracy. The first is the complicated discharge patterns in Lake 
Keowee caused by power plant, hydro and pumped storage operation 
compared to ambient simulation for Cayuga Lake. The second cause is 
that for Cayuga Lake; the assumption of horizontal mixing is more mean- 
ingful since it is a single basin domain, while Lake Keowee is a double 
basin domain connected by a canal. The model, however, can be adapted 
for application to two connected basins. 


SECTION 3 


RECOMMENDATIONS 


The main disadvantage of a one-dimensional thermal model lies in 
the fact that resolution is sacrificed for computational speed. Three di- 
mensional models are bulky and time consuming but have much better 
resolution, however, when long term simulations are necessary, a one- 
dimensional model is recommended. 

The model described here can be modified to include the single 
effects of the various quantities involved in the surface heat transfer 
phenomenon rather than using the equilibrium temperature concept. 

This is particularly recommended for the user who is interested in 
modeling the long term effects of one (for example, evaporation) of 
the quantities involved in the surface heat transfer processes. 

Furthermore, the model can be easily adapted to handle connected 
multiple domains. This recommendation is discussed in the text. 
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SECTION 4 


THE MATHEMATICAL MODEL 


DESCRIPTION 

The one-dimensional variable cross-sectional area model used for 
this study is described in this section. The detailed derivation of the 
governing equations is described in Appendix A, while the boundary 
conditions and the important terms occurring in the governing equations 
are discussed here. 

This model which assumes lateral uniformity was developed to pre- 
dict the vertical temperature profiles for complete annual cycles in a lake. 
The model includes the effects of variation of the horizontal cross-section- 
ai area with depth. Surface heating due to solar radiation absorbed at 
the surface layer and the internal heating due to the transmission of the 
unabsorbed solar radiation to the deeper layers of the lake are taken 
into account. The exchange of mechanical energy between the lake and 
the atmosphere is accounted for through the friction velocity and the 
thermal diffusivity, since these quantities can be related to wind speed. 
Finally, the effects of power plant discharge and intake are considered. 

The full three-dimensional equations of mass and heat balance are: 

|£ = -v^pv (4.1) 

and 

(pC T) = V-pC K-VT - V*pC TV + H (4.2) 

dip p p 

Applying the divergence theorem, and then differentiating and integrat- 
ing, a set of one-dimensional equations are obtained: 

A(z) || = |-^(z)pV^ + lA' (4.3) 

and 

A(z) 1^. (pCpT) =f- 2 [pCpA(z)K^ If] - |-(pCpA(z)TV^) + QA' + A(z)H(z) 

(4.4) 

where, 

z is the vertical coordinate, measured upward from the deepest point of 
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of the lake. As a subscript it marks the vertical component of a 
vector. 

A(*) is the horizontal cross-section of the lake at height z. 
p is the density, 
t is the time. 

is the vertical velocity of flow. 

Cp is the heat capacity. 

T is the temperature. 

K is the vertical component of 'K', the heat diffusivity tensor (including 
^ turbulent diffusivity) . 

H(z) is the source of heat per unit volume. 


A' is ^ 
A IS 


I is the bottom-surface source of mass per unit area. 


Q is the bottom-surface source of heat per unit area. 

The following variables occur in the governing equations and are 
defined thus: 

Density, p, is assumed to vary with temperature in the following form: 

p = 1.02943 + (-2.0 X 10"®)T + (-4.8 x 10"®)T2 gm/cc (4.5) 

Equation (4.5) and the discussion of the constants were given by Sengupta 
and Lick (1974). 

Eddy diffusivity, K , is a function of both thermal and current structure 
of a lake. The form used here was deduced by Rossby and Montgomery 
(1935). 



K 


z 


^ 20^1 


aiRi) 


-1 


(4.6) 


where, 

is the eddy diffusivity under neutral conditions. 

Sundaram et ai. (1969) gave an empirical form for the yearly variation 


K,,, = 0.21 + 0.052 sin (4^ -I- 2.61) cm' /sec* 

ZO aOd 


(4.7) 


6 


A sinusoidal variation of K with a mean 0. 21 cm' /sec and amplitude of 
0.052 cm^ysec is assumed iFPEquation (4.7). 


The semi-empirical constant a. = 0.1 is estimated in this study by 
comparing the values used by Sunaaram et al. (1969) and the original 
value used by Monin and Obukhov (1954). 


The third variable appearing in Equation (4.6) is R|, the Richardson 
number. This is defined as: 



W** 


3T 

az 


(4.8) 


where, 

is the volumetric coefficient of expansion of water, 
is defined as: 

= 0.0 + 1.538 X 10”®(T-4) + (-2.037 x 10”^)(T-4)2 (4.9) 

The constants in the above equation were taken from Sundaram et ai. 
(1970). a can also be estimated from Equation (4.5); however, the 
use of Equation (4.9) leads to a higher accuracy. 

g is the acceleration due to gravity. 

w* is the friction velocity calculated from known conditions above the 
lake. 

The heat source, H, is that part of the solar radiation transmitted 
exponentially through the depth of the lake. The relation used in this 
study has been used by Dake and Harleman (1969) and by Mitry and 
Ozisik (1976) : 

H(z) = 7 i(1-3)A(z) EXP(-n(z-h)) (4.10) 


where, 

8 = fraction of the solar radiation absorbed at the surface, (0.5). 
n = extinction coefficient, (0.75 cm~b. 
h = depth of the lake, (m). 

= net solar radiation reaching the water surface, (cal /cm ^ /day) . 

Equations (4.3) and (4.4) are parabolic and mathematically repre- 
sent a diffusion process. The solution of these equations requires one 
initial condition and two boundary conditions. 
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Initial Conditions 


The temperature of the lake at spring homothermy is taken as the 
initial temperature. 

Boundary Conditions 

Surface: 

where, 

Tg = equilibrium temperature, (®C). 

Tg = surface temperature, (®C). 

Kg = surface heat exchange coefficient, (cal/cm*-s) 

Tg Evaluation 

Brady et ai. (1969) showed empirically that fluctuations in the 
equilibrium, temperature may be conveniently estimated using the approxi- 
mate relationship: 

(4.12) 

where, 

is the dewpoint temperature, estimated from known conditions above 
° the lake. 

Hg is the gross rate of short wave solar radiation. 

Since the dewpoint temperature tends to remain relatively constant through 
a single day. Equation (4.12) indicates that the main source of hourly 
fluctuations in is the solar radiation component. This generally 
reaches a maximuin at solar noon, unless variable cloudiness interferes. 

At nighttime, T. approaches the dewpoint temperature, which acts like 
a relatively invariant datum for periods of 24 hours or less. On an 
annual basis, however, both T . and H are generally much greater in 
summer than in winter. The cmminant ^contribution to the amplitude of 
seasonal fluctuation in T^ is the dewpoint temperature. 

K. Evaluation 
s, 

The form used by Edinger and Geyer (1967) was also used in this 
study: 

Kg = 4.5 + 0.05 Tg + 6f(w) + 0.47f(w) (4.13) 
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3 is found by applying standard curve-fitting techniques to pub- 
lished data pertaining to saturated vapor pressures at various tempera- 
tures; a convenient representation given by Edinger and Ceyer (1967) 
is: 


6 = 0.35 + 0.015 T + 0.0012 T * (mmHgVX) 
m m 


(4.14) 


T + T . 
s d 


m 


(4.15) 


The evaporative wind speed function f(w) used is also similar to that of 
Edinger and Ceyer (1967). 


f(w) = 9.2 + 0.46W2 {VI ^ (4.16) 

where W is the wind speed (m/s) 

Bottom: 

The second boundary condition is at the bottom surface of the lake 
which is assumed to be perfectly insulated. 

lfiz=0 = » 


Numerical Method 

Numerical integration of the governing equations require that these 
equations be replaced by finite difference equations. 

In this study, the following notations are used in the finite differ- 
encin9 of the basic equations: 

Az The thickness of a slice (the basin is divided into equal number 
of slices). This is used as depth increment. 

At The time increment chosen at the n^*^ Integration step. 

n, n-1 The present and last time at which computations were performed. 


J 


The space subscript. This is used as 

T." = T(jAz, t^) = T(z, t) 
where z is the vertical coordinate. 


Equation (4.4) can now be written in difference form as: 
j n _ j n-1 

i 1 \ n-1 1 \ n— 1,.Ai^ n— 1 


At 


n 
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- (-s~) 


1 ,n-l 1 


.n-1 


n-1 


IJ 
T. 
P,A, 

T,l 


jhI(PATV,>i;, - <PATV^)jl,'l 
ir tipAv ,),";,' - ( PAV ,),".-'1 


- [2Tr(h-(jAz)l""’ 


J I 


^ l2.(h-(jAz)l 


n-1 


I P J 


where. 


PjCp J 


(pAK^AT)";'^’ = pAKJ(j+i)Az, ‘ tP"^1 

(pAK^AT)jP^ - pAK^[(j-i)Az, ViHT""^ - tP’^I 


9*i _1 gT' 

K = K [1 + a,(- —A) ' (a ).z. ^ 
z. - zo 1' )hr^' v^l J 9z 


9T1 


1+1 


= K^t(H)4lH1 + a,(- )~^-(»y)j[}(TP'' +tJ’-’)KH 4 i: 


*Vj ' A, + AjC. - «) + Aj(T. - 4)> 


n-1 


-n-1. 


P,.j = 3, * Tjlj- . = 


A =J!^I 
n-J h ' 




Terms like AT.., and AT. , can be defined as 
1+1 j-i 
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The sequence in which calculations are performed is as follows: 

The dependent variables, T, K , w*, a , p, Tp, and K , are 
initialized. The area of each slice il calculared. Tn& title o^the pre- 
sent year Is printed. The values of the variables, K , w*, a , p, Tp, 
and K , are then calculated. The temperatures of the slices are final^ 
calcul^ed. If the temperature profile is unstable, n^ixing of the unstable 
portion of the profile is undertaken. ' 

During the next time step, the temperatures are updated, and the 
dependent variables are calculated again. 

The values of the temperature, T, eddy diffusivity, K^, number of 
days and the surface heat transfer coefficient, K^, are printed every 
time step, every day or normally at the end of each month. At the end 
of the present year, the new-year title is printed and computations con- 
tinue as listed above. A detailed description of the programs and flow 
chart has been published in a user's manual, Lee and Sengupta (1980) . 
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SECTION 5 


MODEL CALIBRATION 


The model described in the previous section was appiied to Cayuga 
Lake, New York. The main reason was to calibrate the model. Cayuga 
Lake was chosen for this calibration because of the vast amount of nu- 
merical work done on this lake. These include Sundaram and Rehm (1971), 
Sundaram et al. (1971), Sundaram and Rehm (1972), and Mitry and Ozisik 
(1976). The model was also verified using Lake Keowee data base; this 
is described in the next chapter. 

CAYUGA LAKE APPLICATION 

Cayuga Lake is one of the finger lakes in west central New York. 

It is the longest, about 40 miles (64 km), of the group, and the second 
largest, about 66 square miles (171 sq. km), with steep banks cut by 
inflowing streams. It averages 2 miles (3 km) wide about 500 ft (150 
meters) deep, with a surface elevation of 381 ft (115 meters). Tang- 
hannock Creek and its 210-foot (64-meter) cataracts flow into the lake. 

The lake is drained to the north by the Seneca River, which joins the 
Oswego River and flows into Lake Ontario. Cayuga Lake is connected 
by canal with Seneca Lake on the west and Erie Canal on the north. 

The city of Ithaca lies at the’ southern end of the lake. 

INPUT QUANTITIES 

In the numerical integration of the governing equations for Cayuga 
Lake, the following relations were used. The numerical values were 
taken from Sundaram et al. (1971) and Mitry and Ozisik (1976). 

Density, p, varied as in Equation (4.5), Section 4, of this report. 

Eddy diffusivity varied as in Equations (4.6) and (4.7) with 
=0.1 and the Richardson number as given in Equation (4.8). 

The volumetric coefficient of expansion of water, a , used is shown 
in Equation (4.9). ^ 

A sinusoidal relation was assumed for the variation of the friction 
veiocity, w*, with time; normally, w* is defined as 


w* = -ytTg/p) 


(5.1) 


The form used is 


w* = sln(|g + C^) (5.2) 

where, 

= Average value of w*, (3.048 cm /sec). 

= Half the annual variation of w*, (0.762 cm /sec). 

C- = Phase angle, (2. 61) (chosen In such a way that at time t = 0, 
w* = initial value of the friction velocity). 

The heat source, -H, is taken to vary as in Equation (4.10), with 
3 = 0. 5, n = 0.75 cm and h = 70.0 m ( = 200 ft). 

The net solar radiation, <l> , reaching the water surface is taken 
from Mitry and Ozisik (1976): ° 

♦„ = Aj + S2Sin(^t+Cj) (5.3) 

where, 

* 

— 3 

A^ = Average value of t(>^, (6. 14 x 10 cal/cm^s). 

~ 3 

B^ = Half the annual variation of (3.52 x 10 cal/cm^s). 

C 2 = Phase angle, (0.049) (chosen in the same way as C.j). 

The heat flux, Q, in Equation (4.4) is defined as 

Q = (pCpATQp)/(A(2)) (5.4) 

where. 


Qp = Volumetric discharge from the power plant. In this study 

Qp - 1.508 X 10° cm^/sec; this value is chosen to correspond 
ter the pumping velocity used by Sundaram et al. (1971), (1/4 
ft /day) . 


AT = 10°C = Assumed temperature change through the condensers 
of the power plant. 

The pumping velocity is defined as 


V 


A(z) 


(5. 5) 


The pumping velocity term effects are felt between the intake and the 
level at which the heated effluent becomes neutrally buoyant (effective 
discharge level). 
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- 4 

K = Assumed constant at 5.65 x 10 cal/cm^-s. 
s 

-4 

For a postulated 3500 MW plant for Cayuga Lake, a 8. 79 x 10 
cal/cm^s of waste heat will have to be rejected. 

Surface area of Cayuga Lake is 66 sq. miles (171 sq. km). 

The intake to the power plant is fixed at 38.1 m 125 ft) from 
the surface of the lake. 

Initial temperature is taken as the temperature of the lake at spring 
homothermy (usually coincides with the minimum surface temperature 
attained by the lake). For Cayuga Lake, this occurs around March, 
and the temperature assumed is = 2. 9®C. 

Finally, a sinusoidal variation is assumed for the variation of the 
equilibrium temperature with time. 

Tg ^Aj + BjSinC^t+Cj) (5.6) 

where, 

Aj = Average value of T^, (11®C). 

Bj = Half the annual variation of T^, (16®C). 

C^ = Phase angle, chosen in the same way as C.| or C^* 

The justification for using Equation (5.6) is described in Appendix 
D of this report. 

Two topographies are studied in the application of this model to 
Cayuga Lake; these are a right circular cylinder and a circular para- 
boloid approximations for the lake. 

Cylindrical Topography 

5 

The area (radius, B = 7. 38 x 10 cmL^f the lake is assumed con- 
stant throughout the depth. Term A' or = 0; see Equations (4.3) 
and (4.4). 

Circular Paraboloid Topography 


The lakeJs assumed to be a circular paraboloid with surface radius, 
B - 7. 38 X 10^ cm. The area at any depth, z (measured from the deep- 
est point of the lake) is given by: 



(5.7) 


Thus, A' is a constant: 
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RESULTS 


A' 

^ " h 


(5.8) 


Computations for a yearly cycle for Cayuga Lake are presented. 

The verification data base consists of vertical temperature profiles com- 
piled by Henson et al. (1961). The comparison of simulated and observed 
vertical temperature profiles are shown in Figures 1, 2 and 3. Each 
figure shows five profiles representing observed (M), discharge (cylin- 
drical domain - CD) ; discharge (parabolic domain - PD) ; no-discharge 
(cylindrical domain - CN); and no-discharge (parabolic domain - PN). 

The no-discharge simulations are in good agreement with the data 
(the data was for no-discharge conditions) . 

The parabolic case has somewhat better agreement since it represents, 
qualitatively, the decrease in area with depth. However, the closeness 
of the simulated results for the two cases is surprising. Most lakes have 
the rate of decrease in area with depth greater than a paraboloid, which 
has a linear decrease. Thus, when realistic area changes are used, a 
greater difference between cylindrical and paraboloid cases can be ex- 
pected. 

The effect of discharge is significant only in the top layers till July. 
This is because the heated discharge rises to the surface during the first 
half of the year. For the later months, the discharge temperature is 
lower than the surface temperature causing the discharge to reach static 
equilibrium somewhere below the surface. Thus, significant thermal 
effects of discharge are seen at mid-depths till December. The tempera- 
tures are higher at these depths for the paraboloid topography. In 
general, a temperature difference of the order of 3®C over the no-dis- 
charge case, can be seen. At the end of the annual cycle, a residual 
temperature Increase of 1.75®C is detectable. 

In the above paragraph, the discharge from the power plant is 
treated as a plane source which is injected into the domain at a level 
where the discharge temperature equals the temperature of that level. 

This level could be called the effective discharge level. The effects of 
the pumping velocity term, therefore, can only be felt from the intake 
level to the effective discharge level. A temperature rise of 10®C through 
the condensers was assumed. This temperature difference justifies the 
use of density as a function of temperature. 

Figures 4 and 5 show the annual stratification cycle. From these 
figures it is evident that the surface temperature difference between the 
four cases is less than 2°C over the annual cycle. However, at mid- 
depths the paraboloid discharge case shows a 5°C difference compared 
to ihe no-discharge case; for the cylindrical discharge case this difference 
is 3X. Generally, the highest surface temperatures are reached after 
150 days, while the highest equilibrium temperature occurs after 120 days. 
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Thus, there is approximately a 30 day lag in the surface temperature 
response. The maximum temperatures at mid-depth occur after 240 days 
for the no-discharge case; for the discharge case, the corresponding time 
is 210 days. No significant phase lag between cylindrical and paraboloid 
cases are observed. 

Figures 6 and 7 show the eddy diffusivity variation with depth 
(time in days used as parameter) for the cylindrical and paraboloid 
cases respectively. It is observed that thermal discharge causes in- 
crease in eddy diffusivity in the epilimnion owing to increased mixing. 

No significant changes are seen in the hypolimnion. There is a temporal 
increase in the differences between the discharge and no-discharge cases. 
These observations are true for both domains. However, the diffusivity 
values are larger for the paraboloid cases. Also, at any given time, 
the paraboloid case shows deeper thermoclines. 
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Figure 1. Vertical temperature profiles, Cayuga Lake 
(from 0 to 90 days) 
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Figure 2. Vertical temperature profiles, Cayuga Lake 
(from 120 to 210 days) 
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Figure 3. Vertical temperature profiles, Cayuga Lake 
(from 240 to 330 days) 
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Figure 4. Stratification cycle cylindrical domain, Cayuga Lake 
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Figure- 5. Stratification cycle parabolic domain, Cayuga Lake 
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Figure 7. Variation of eddy diffusivity with depth paraboloid domain, 
Cayuga Lake 
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SECTION 6 


MODEL VERIFICATION 


The model calibrated in the last section was applied to Lake Keowee. 
There were two objectives. The first was to simulate the stratification 
behavior of Lake Keowee for the period 1971 to 1979, the entire existence 
cif this river-dammed lake, which started receiving power plant heated 
discharges in 1973. The second objective was to test the accuracy of 
this model for a situation where horizontal uniformity (see Appendix A) 

0f temperature was somewhat suspected since there are two basins con- 
nected by a canal. 

LAKE KEOWEE APPLICATION 

Lake Keowee is located 40 km west of Greenville, South Carolina. 

It is the source of cooling water for Oconee Nuclear Power Station (ONS). 
It was formed from 1968 through 1971 by damming Little and Keowee 
Rivers. A connecting canal (maximum depth 30. 5 m) joins the two main 
arms of the lake. Flow out of the lake, except for negligible leakage 
through Little River dam, is through the Keowee Hydro Station. Lake 
Keowee also exchanges water with Lake Jocassee-pumped storage station. 
The three-unit steam electric ONS with a net capacity of 2580 Mwe 
started operating in July 1973. ONS operated at annual gross thermal 
capacity factors (GTCF) of 11, 28, 69 and 59% in the years 1973 through 
1976, respectively, Duke Power Company ( 1 976) . From 1977 through 
1979 the GTCF varied from 65 to 75 , Duke Power Company (1979). 

INPUT QUANTITIES 

Since long term effects are the important factors in this sutdy, the 
input data to the one-dimensional model were mainly monthly averages 
supplied by Duke Power Company (1976 and 1979). These are: 

1. Maximum depth of Lake Keowee = 45 m (si 50 ft). 

2. Initial temperature (minimum surface temperature of Lake Keowee) 

= 7.5®C. 

3i ONS condenser cooling water flowrate, Q , Table 1. 

4. Temperature rise in condenser cooling water. Table 2. 

5. Outflow through Keowee Hydro Station, Table 3. 
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6. Flow through Jocassee-pumpted storage station. Table 4. 

7. Wind speed. Table 5 (0.2 meters from surface). 

8. Solar radiation reaching the water surface, <j)^. Table 6. 

9. Dewpoint temperature. Table 7. 

RESULTS 

Computations were carried out for nine years, 1971 through 1979, 
using the above monthly-averaged data. The measured temperature 
profiles referred to in this discussion are horizontal -averaged data col- 
lected from seven different locations in Lake Keowee, as shown in Figure 
8. This is in approximate agreement with the lateral uniformity assumed 
in the formulation of this and other one-dimensional models. 

Computed (with discharge) and measured temperature profiles for 
1971 to 1979 are presented together whenever possible. Figures 9 through 
17. A second set of figures, 18 to 26, for the nine simulated years show 
the computed temperature profiles (with and without discharge). Other 
figures presented include: a typical variation of eddy diffusivity with 

depth (time in days as parameter). Figure 27; and a typical temporal 
variation of the depth of the thermocline. Figure 28. Finally, the 
variation of the surface heat exchange coefficient with time in years 
(Figure 29), and the stratification cycle for Lake Keowee (1971-1979), 
Figure 30, are presented. All the above figures are computer drawn. 

For 1971, the year the lake came into existence. Figure 9 shows 
the computed temperature profiles (solid lines) and the measured tem- 
perature profiles (broken lines). The surface temperature comparisons 
between the measured and computed are excellent except for the start 
of the heating season when the computed temperatures seem to lag the 
measured. The maximum lag occurs in April and reduces with the 
heating season and completely disappears in October when cooling is in 
full swing. The hypolimnetic temperatures show good agreement with 
measured data. The mid-layer temperatures compare perfectly for the 
first half of the year. From August to November, the trend changes 
and the computed profiles lead the measured. However, the modeling 
of the formation, maintenance and decay of the thermocline, is good. 

This can be seen by comparing the shapes of the measured and computed 
temperature profiles. The comparison improves with increasing years. 
Figures 10 to 17, and generally follows the same trend as outlined above. 
The comparison of measured and computed temperature profiles cannot 
be expected to be more precise for the following reason. The measured 
values were taken on some specific day (different for some months) 
during the month while the computed results show the temperature pro- 
files at the end of the month. 

To estimate the impact of ONS operations on the thermal dynamics 
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of Lake Keowee, a set of simulations (with and without discharge) are 
presented. These are summarized in Figures 18 to 26. Two years are 
selected for detailed discussion, 1975 and 1978. The earlier year. 

Figure 22, was a year of very high ONS annual gross thermal capacity 
factor and the later year. Figure 25, is representative of the situation 
towards the end of computations. In Figure 22, although the discharge 
surfaced only during the early and later months of the year, the effect 
of discharge (broken lines) is clearly evident during most of the months 
at the subsurface layers (especially below the thermcline) . Figure 25 
shows the same situation for 1 978. This time, however, the discharge 
surfaced from January through May and then from October through 
December, at total of three months more than in 1975. The hypolimnetic 
temperatures were ciso higher. Finally, by the end of computations, 
1979, there was a residual temperature difference of aimost 3°C due to 
discharge. ONS started operations in July 1973, and accordingly, the 
first three years. Figures 18 to 20, show no differences between dis- 
charge and no-discharge simulations. 

Figure 27 shows the variation of the eddy diffusivity with depth 
for 1975, one of the nine years simulated. In this figure, the maximum 
diffusivity always occurs at the surface for every month. The non- 
stratified months show a constant diffusivity while the stratified months 
show very high mixing in the epilimnion and very low mixing in the 
hypolimnion. 

Figure 28 shows the temporal variation of the depth of the thermo- 
cllne for 1975. This figure Is In agreement with one of the, generally 
accepted theories of the formation, maintenance and decay of the thermo- 
cline. 


Finally, a summary of the stratification cycle for Lake Keowee for 
the nine years simulated is presented in Figure 30, and the variation of 
the surface exchange coefficient for the same period is shown In Figures 
29 and 30. The maximum surface temperatures occur around 240 days, 
and the minimum surface temperatures around 60 days for each year. 

The equilibrium temperature of the lake is superimposed. The maximum 
equilibrium temperatures occur around 210 days. The surface temperature 
therefore lags the equilibrium temperature by 30 days. Two curves are 
Shown for the mid-layer temperatures to highlight the effect of ONS 
operations. The temporal increase of this effect can be seen clearly. 

The maximum difference occurs during the last quarter of each year. 

The variation of the surface heat exchange coefficient, as expected, 
follows the same trend as the variation of the surface temperatures. 


26 


TABLE 1. Oconee Nuclear Station Condenser Cooling Water Flowrate (m^/min) 


Month 

1973 

1974 

1975 

1976 

1977 

1978 

1 979 

January 


3069. 3 

4612.4 

6069. 3 

5045. 8 

6176.7 

7207. 7 

February 


3069. 4 

3694.9 

4440.2 

4985. 2 

6444. 6 

731 9.9 

March 


2976. 9 

5456. 8 

4874. 3 

51 1 3. 5 

51 95. 7 

741 9. 5 

April 

— 

2807. 3 

5570. 8 

4272.1 

601 3. 6 

4811.8 

7275. 8 

May 

— 

2164.6 

6494. 3 

3970. 7 

6302. 4 

4984. 2 

41 89. 1 

June 

— 

4171. 8 

6574. 2 

51 97. 6 

4385. 3 

5659. 9 

5381.2 

July 

1 890. 2 

5334.6 

7104.2 

5830.0 

5038. 6 

7058.8 

4753. 3 

August 

1910. 3 

i 

4727. 1 

7510.1 

7248. 3 

5708. 9 

791 4.9 


September | 

2099. 5 

5961.4 

7201.6 

6785. 4 

6964.0 

6557. 3 

— 

October 

2232.5 

4953. 4 

6993.4 

5637. 8 

6754. 7 

7407. 4 


November 

2170.7 

4202. 1 

7467. 1 

5809. 2 

4697. 6 

6065. 1 

— 

December 

32 84. 6 

5225.6 

6850. 9 

491 4. 8 

5854. 6 

6503. 5 

— 












TABLE 2. Oconee Nuclear Station Condenser Temperature Rise, AT(**C) 


Month 

1973 

1974 

1975 

. 1976 

1977 

1978 

1979 

January 

— 

mm 

6. 3 

Dl 

12.5 

9.0 

10.3 

February 

— 

IB 

4.8 

B 

11.4 

11.0 

10.4 

March 

— 

8.4 

8.2 

7.1 

10.4 

13.2 

9.6 

April 


8.0 

8.3 

5.1 

11.4 

9.7 

9.9 

May 

— 

B 

8.8 

5.8 

9.4 

10.1 

8.2 

June 



8.8 


8.4 

8.1 

7.1 

July 

5.3 


8.3 


7.4 

7.9 

5.0 

August 

4.6 

4.8 

7. 8 

8.5 

5.0 

7.5 

— 

September 

5.3 

5.8 

B 

8.0 

5.0 

7.6 

— 

October 

7. 3 

3.5 


7.8 

3.8 

6.2 

— 

1 

November 

7.7 

7.9 


6.7 

6.2 

1 

8.4 

1 

— 

December 

4.1 

5.9 


8.4 

7. 9 

1 

7.2 













TABLE 3. Monthly Average Flowrates (m^/sec) - Lake Keowee Hydro Station 


Month 

1971 

1972 

1973 

1974 

1975 

1976 

1977 

1978 

January 

— 

17. 5 

22.2 

71.2 

27.0 

39. 5 

35.5 

30.1 

February 


32.9 

17. 2 

118.1 

50.8 

36.2 

33.6 

22.6 

March 

— 

9.0 

33.6 

20.7 

82.7 

1 

42.4 

25.4 

25.7 

April 

— 

3.0 

30.1 

16.7 

42.5 

40.7 

50.5 

18.7 

May 

— 

21.2 

31.6 

35.9 

63.8 

53.3 

21.6 

52.3 

June 

— 

22.9 

55.7 

11.8 

16.1 

53.2 

27. 5 

43.1 

July 


14.3 

9.7 

43.1 

20. 3 

33.2 

20.1 

12.2 

August 

— 

15.8 

16.1 

28. 8 

39. 9 

10.8 

55.2 

28.7 

September 

22. 6 

20.2 

18.4 

17.0 

35.1 

16.0 

16.5 

33.1 

October 

1 4. 9 

6.1 

3.3 

5.1 

49.8 

39.2 

7.5 

25.2 

November 

18.5 

30.8 

5.6 

23.0 

47.2 

27.2 

22.5 

13.4 

December 

8. 8 

17. 3 

19. 3 

12.4 

44.1 

30.1 

20.5 

13.4 

.Annual Average 

16.2* 

17.6 

21.9 

H 

43. 3 . 

ID 

27.2 

32.1 


* Average of S;iptember through December. 
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TABLE 4. Lake Jocassee Hydro FFows (cfs) 


Month 

1973 

1974 

1975 

1976 

1977 

1978 

1979 

January 

— 

: ; 

— 

505 

-608 

187 

1 

February 

— 

— — 

1 : 

509 

206 

54 

— 

March 

— 

— 

j 

562 

1181 

516 

— 

April 


— 


-79 

1046 

370 

— 

May 

• 

— 

— 

555 

27 

1762.2 

— 

June 

— 

— 

— 

1032 

489 

180.2 

— 

July 

— 

— 

— 

112 

-831 

— 

— 

August 



-163 

-747 

-460 

— 


September 

— 

— 

554 

1 

39 

101 

— 


October 

— 

— 

765 

304 

152 

— 

— 

November 

— 

i 

615 

-251 

448 

1 


December 

— 

— 

301 

-915 

461 

1 

— 
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TABLE 5. Lake Keowee, Wind Speed (cm /sec) 


Month 

1971 

1972 

1973 

1974 

1975 

1976 

1977 

1978 

1979 

January 

6. 69 

6. 69 

7. 22 

5.8 

6. 3 

7.4 

8.04 

7.9 

8.6 

February 

9.3 

9.26 

7.3 

5.8 

7.6 

8.5 

8. 4 

6.8 

7.2 

t 

March 

9.2 

9. 20 

7.1 

7.7 

9.6 

7. 9 

7.7 

7.6 

7.9 

April 

8.72 

8.7 

8. 44 

8.7 

7.6 

7.6 

7.6 

7.6 

7.6 

May 

7.5 

7. 53 

6.83 

6.8 

4.8 

7.3 

6.2 

6.7 

6.7 

June 

5. 65 

7. 95 

3.04 

6. 96 

5. 82 

6. 4 

6.7 

4.7 

4.8 

July 

6. 48 

6.64 

5. 32 

5.2 

5.1 

5.9 

5.8 

5.7 

— 

August 

5. 75 

6.07 

5.1 

5.87 

5.4 

6.7 

5.4 

5.1 

— 

September 

5. 77 

5. 47 

6.8 

6.74 

7.3 

7.13 

5.3 

5.7 

— 

October 

7.02 

7.17 

7.1 

5.7 

7.7 

7. 21 

7.2 

6.6 

— 

November 

7. 53 

7.13 

8.14 

CM 

• 

6.9 

7. 27 

7. 5 

5.8 

— 

December 

8.3 

6.8 


6.9 

7.2 

8.2 

7. 2 

7.3 

1 

— 

















TABLE 6. Lake Keowee Gross Solar Radiation (Langleys) 


Month 

1971 

1972 

1973 

1974 

1975 

1976 


1978 

1979 

January 

167.0 

176.0 

162.7 

191.4 

191.4 

209.8 

1977 j 
205.5 

227.0 

208.0 

February 

264.4 

257.6 

279.5 

226.9 

226.9 

310.9 

317.6 

308.0 

251.0 

March 

264.4 

352.5 

348.5 

326.1 

326.1 

338.6 

328.5 

408.0 

373.0 

April 

457.0 

448.0 

449.3 

397.7 

397.7 

496. 9 

427.3 

429.0 

479.0 

May 

480.5 

433.6 

449.5 

436.0 

436.0 

448. 4 

473.0 

513.0 

513.0 

June 

478.0 

564.3 

507.7 

559.3 

559.3 

480.2 

543.3 

598.0 

— 

July 

409.0 

493. 8 

496.9 

459.5 

459.5 

288. 3 

551. 8 

568.0 

— 

August 

428.2 

453.5 

391.6 

480.0 

480.0 

480. 4 

423.9 

461.0 

— 

September 

329.0 

386. 3 

338.4 

339.2 

339. 2 

345.1 

350.7 

385.0 


October 

261.3 

298.1 

341.7 


302.5 

287.5 

286. 6 

369.0 

— 

November 

247.7 

220.9 

247.6 

231.1 

231.1 

237.5 

1 96. 2 

232.0 

— 

December 

147.7 

148.0 

154.0 

181.9 

181.9 

195.0 

178.2 

191.0 
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TABLE 7. Lake Keowee Dewpoint Temperature (*’C) 


Month 

1971 

1972 

1973 

1974 

1975 

1976 

1977 

1978 

1979 

January 

3.0 

1 

1.67 

1.0 

8. 2 

1 

3.0 

-1.0 

-6.6 

-2.8 

-3.33 

February 

0.9 

-2.22 

1 

-1.0 

0.0 

3.5 


-2.78 

-5.0 

0.0 

March 

6. 3 

1.11 

10.0 

6. 3 

2.2 


6.0 

1.2 

5.0 

April 

7.5 

6.6 

7.7 

10.7 

7.2 

11.2 

10.2 

9.6 

9.2 

May 

17.2 

11.11 

14.3 

17.2 

17. 5 

14.0 

15.4 

14.0 

14.0 

June 

18. 8 

13.13 

20. 25 

17.8 

19.0 

18.3 

18.0 

19.4 

19. 4 

July 

20.0 

18.77 

22. 2 

21.0 

21.3 

19.8 

20.2 

20.8 

20.7 

August 

19. 44 

22. 22 

21.7 

21.0 

21.0 

18.0 

20.7 

20.8 

21.0 

September 

18. 33 

18.8 

20.8 

17. 5 

16. 2 

15. 4 

18.7 

15.5 

15.7 

October 

13.88 

11.5 

13. 5 

10. 2 

12.4 

8. 2 

9. 2 

9. 3 

10.0 

November 

2.88 

5.9 

7.:2 

6.0 


1.0 

7.0 

9.0 

10.3 

December 

5.5 

4. 0 

3.2 

3.8 


-1.5 

0.4 

0.4 

1 

0.9 
















Figure 8. Map of Lake Keowee 
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Flqus^e 9. Measured (broken lines) and predicted temperature profiles 
Lake Keowee, 1971 
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Figure 10. Measured (broken lines) and predicted temperature profiles 
Lake Keowee/ 1972 
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Figure 11. Measured (broken lines) and predicted temperature profiles 
- Lake Keowee, 1973 
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Figure 12. Measured (broken lines) and predicted temperature profiles 
Lake Keowee, 1974 




Figure 13. Measured (broken lines) and predicted temperature profiles 
Lake Kecwee, 1975 
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Figure 14. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1976 
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Figure 15. Measured (broken lines) and predicted temperature profiles. 
Lake Keowee, 1977 






OO’o'l 

, 01 * 


00 



tiJ JHid3Q 








} 

t 

'I 


42 


Figure 16. Measured (broken lines) and predicted temperature profiles 
Lake Keowee, 1978 
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Figure 17. Predicted discharge temperature profiles. Lake Keowee, 1979 
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Fiqure 18. PrecRcted hQ*“discharae temperature profiles^ Lake Keowee, 1971 










Figure 19. Predicted no-discharge temperature profiles. Lake Keowee , 1 972 
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Figure 20. Predicted discharge (horken lines) and no-discharae temoerature oroflles 
Lake Keowee, 1973 
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Figure 22. Predicted discharge (borken lines) and no-discharoe temoerature orofiies 
Lake Keowee, 1975 
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Figure 23. Predicted discharge (borken lines) and no-discharge temperature profiles 
Lake Keowee, 1976 
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Figure 2-4. Predicted discharge (broken lines) and no-discharqe temperature profiles 
Lake Keowee, 1977 ' ' 
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Figure 26. .Predicted discharge (broken lines) and no-discharge temperature profiles 
Lake Keowee, 1979 
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Figure 27. Variation of eddy diffusivity with depth. Lake Keowee, 1975 
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Figure 28. Monthly variation of the depth of the thermocline. 
Lake Keowee, 1975 
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Figure 29. Temporal variation of the surface exchange coefficient. Lake Keowee, 1971-1979 
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Figure 30. Stratification cycle. Lake Keowee, 1971-1979 
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APPENDIX A 


DERIVATION OF THE MATHEMATICAL MODEL 

DERIVATION OF THE VARIABLE CROSS-SECTIONAL AREA MODEL 

Assuming lateral uniformity (horizontally uniform lake) , the isotherms 
coincide with the isopycnals. The equations for such a mt^ei are obtained 
from the basic balance equations of mass and heat: 

If = -V • pv (A.1) 

|t(pC„T) = 7 • pC R • V • pC TV + H (A, 2) 

tfc p p p 

where, 
p = density, 
t = time. 

V = velocity of flow. 

Cp = heat capacity. 

T = temperature. 

K = heat diffusivity tensor (including turbulent diffusivity) . 

H = source of heat per unit volume. 

Most models have started with the above equations and neglected 
all the terms involving the two horizontal components (or dimensions). 

The negiect of the horizontal divergence in such simple models is not 
always justified. There are at least two reasons for the existence of 
horizontal divergence in real lakes. 

1. The variation of horizontal cross-sectional area of the lake with depth. 

2. The existence of sources of heat and matter efflux on the bottom of 
the lake at depths above the deepest point. 
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The need to include these in the diffusion equations of lakes was 
already felt; Lerman and Stiller (1969), Dutton and Bryson (1962) and 
Tzur (1973). Only Tzur (1973) formulated corrected diffusion equations. 

Another complication in applying a 1-D model to a lake is caused by 
seiches and internal waves. When these are present, the height of an 
isopycnal loses its meaning. According to Sweers (1968), averaging 
I simultaneous temperature profiles at different points of a lake or taking 

I profiles during a short interval of time in one or more points can lead 

I to distortion of the shape of the thermocline. He recommends averaging 

i the depths of isotherms instead. In a lake of constant depth this is 

probably the best available approximation but fails for a lake of varying 

depth. 

^ The effects of area change with depth are included by the following 

1 treatments of Equations (A.1) and (A. 2) . Integrating Equation (A. 1) 

over the volume of water below height h measured from the deepest point 
in the lake yields 

J dV = / V • pVdV 

. V 

; Using Gauss theorem on the right hand side yields 

/f£dV = -p/f^ • VdS (A. 3) 

i S 

where S is a surface completely surrounding the volume V, hence 

dS = dC + dA 

where C is the surface area of the part of the bottom of the lake that 
is bounded by the contour at height z. As a subscript, z marks the 
vertical component of a vector. 

Using dV = Adz in Equation (A. 3) yields 

h«^ 

/-It Adz =-p/fi • VdC - p/f1 • VdA 
o^* C A 

I • 6 • ' . 

'?|fAdr=-?pV„dC-pV/(h) 
o o 

Integrating Equation (A. 2) over the volume of water below height h 
measured from the deepest point in the lake yei Ids 
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J*l*pC^TdV = J-(V -pC^R'VTjdV - / V-pC TVdV + / HdV 
y«iP y P V ” V 


Applying the divergence theorem to the first two terms on the right 
yields 


M(z)|^pCpTdr = ff\- (pGpR. 7T)dS - /(ft-pCpT\?)dS 

h 

+ I A(z)H(z)dz 
o " 

Using dS = dA + dC yields 


I' A(z)|;-pC Tdz = / (pc K- V T) dA + / (pC K- 7 T) dC 
0(Jip ©P qP 

h _ h 

- f (pCpTV)^dA - / (pCpT^)^dC 


i.e. 


h 

+ / A(z)H(z)dz 
o 


/A(z)|;rpC Tdz = pC A(h)IK- VT], + f p C [K- 7T1 dC 

jj«p p z^cc n 


- pCpA(h)TV^ - / P^CJ^V^dC 


h 

+ /A(z)H(z)dz (A. 5) 

o 

where, 

Z the vertical coordinate measured upward from the deepest point of 

the lake. As a subscript, it marks the vertical component of a vector. 

C = the surface area of the part of the bottom of the lake that is bounded 
by the contour at height z. 
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n = subscript that marks the component of a vector that is perpendicular 
to the iake-bottom positive upwards. 

A(z) = the horizontai cross-section of the lake at height z. 

As ail the properties of the lake are constant over the horizontal 
isopycnals, the integration over the isopycnals becomes simple multipli- 
cation by their areas. 

Differentiating Equations (A. 4) and (A. 5) with respect to the height, 
a set of 1-0 equations are obtained from Equation (A. 4): 




3C 
nd z 


(A. 6) 




The last term in the above equation is the mass addition term (source 
term) from Equation (A. 5): 




. lp^C^(R . 7T )„ - P^c . AH (A. 7) 

in the above equation, 

p C (K • VT) = Conduction 
c c n 

Heat addition terms 

p C T V = Convection 

c c c n 

Because the horizontal gradients vanish. Equation (A. 7) can be simplified 
further by noting that 

(R-VT)^=K^U (A. 8) 

Following a similar treatment used by Y. Tzur (1973), the slope of the 
bottom of lakes is usually small; it rarely exceeds a gradient of 0.1, 
Assuming this and considering an elemental area shown below, 

dAfl 
dZ 
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Un e = = gradient of the lake bottom 

Cos 6 = Cos [tan” ^ (A. 9) 

also 

Cos 0 — (A. 10) 

From Equation (A. 9) and (A. 10) 

^ sCosCtan"^ (A. 11) 

- Cos[arctan (gradient)] 

ST) - average of Cos arctan 
gradient of the bottom surface 
at that depth. 

From Equation (A. 11) 

dA = ndC 


or 


dA _ 

ar - 



i.e. 

dC _ 1 dA 
dz “ n dz 

Since Cos arctan (0.1) = 0.995, n s — = l. 

ri 

Aiternatively, n can be incorporated into the surface sources (equation 
(A. 7)). 


The 1-D equations can now be written in more concise and conveni- 
ent forms: 

A|f = -fjApVj, + lA' (A.12) 

and 


A|f(pCpT) 


* QA' + AH 


(A. 13) 


where, 
A' 


dA 

dz 
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I — bottom-surface source of mass per unit area. 

Q = bottom-surface source of heat per unit area. 

Boundary Conti tions 
Surface: 

Equilibrium Temperature Concept 

The bulk temperature of a vertically-mixed body of water under 
natural conditions, Tn, tends to increase or decrease with time, accord- 
ing to whether the sum, SH, of its heat inputs (net solar and atmospheric 
radiation, thermal discharge and heat outputs; back- radiation, evaporation, 
and condensation) is positive or negative: 

= EH/pCph (A. 14) 

Following a similar procedure used by Edinger and Ceyer (1967), Equa- 
tion (A. 14) can be transformed (without linearizing the temperature 
dependence of the components of EH) to yield 


"ar " ^ n^ ^P^p^ ^ 

where, 

K_ = surface exchange coefficient that depends on the water temperature 
” and wind speed 

T- = equilibrium temperature, defined as the hypothetical water surface 
^ temperature at which the net rate of surface heat exchange would 
be zero. 



together via water temperature and meteorological conditions. For this 
reason iterative techniques are often used when applying Equation (A. 15) 
in a predictive roile and because the rate of convergence is extremely 
rapid. The first step often yieids satisfactory results. 

Brady et al. (1969) have shown empirically that fluctuations in the 
equilibrium temperature may be conveniently estimated using the approxi- 
mate relationship: 


T,. = T. 


H,/K^ 
s s 


(A. 16) 
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where, 

Td = dewpoint temperature. 

H = gross rate of shortwave solar radiation. 

9 

Since the dewpoint temperature tends to remain relatively constant 
through a single day. Equation (A. 16) indicates that the main source of 
hourly fluctuations in T- is the solar radiation components. This gene- 
rally reaches a maximum^ at solar noon, unless variable cloudiness inter- 
feres. At nighttime T_ approaches the dewpoint temperature which acts 
like a relatively invariant datum for periods of 24 hours or less. On 
an annual basis, however, both T . H are generally much greater in 
summer than in winter. The dominanf contribution to the amplitude of 
seasonal fluctuation in T^ is the dewpoint temperature. 

Exchange Coefficient Evaluation 

The same form used by Edinger et al. (1967) was used: 

K_ = 4.5 + 0.05 T_ + 6f(w) + 0.47 f(w) (A. 17) 

where. 


T. = surface temperature. 

3 is found by applying standard curve-fitting techniques to published 
data pertaining to saturated vapor pressures at various temperatures; 
a convenient representation given by Edinger and Geyer (1967) is 


where. 


3 = 0.35 + 0.015 + 0. 0012 T^^^^ (mmHg/C) (A. 18) 





(A. 19) 


The evaporative windspeed function f(w) used is also similar to that 
of Edinger and Geyer (1967): 

f(w) = 9.2 + 0. 46 W*(Wm"^mmHg‘^ (A. 20) 


where, 

W = wind speed (m/S) . 

The first boundary condition is at the surface and can now be stated as: 

‘•s = '>S^tr|z=h = Ks'Te - V 
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Bottom : 


The steond boundary condition is at the bottom of the iake which is 
assumed to be perfectly insulated 


3T 

3z |z=0 


= 0 


(A. 22) 


and the initial condition is the temperature of the lake at spring homo- 
thermy. 


^initial ” "^o 


(A. 23) 
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APPENDIX B 

MEASURED TEMPERATURE PROFILES 
(Individual Stations) 
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Figure B-1. Lake Keowee measured temperature profiles, 1971 - Station 500 
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Figure B-2. Lake Keowee measured temperature profiles, 1972 - Station 500 
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Figure B-3. Lake Keowee measured temperature profiles, 1973 - Station 500 
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Figure B-4. Lake Keowee measured temperature profiles, 1974 - Station 500 
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Figure B- 5. Lake Keowee measured temperature profiles, 1975 - Station 500 
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Figure B-6. Lake Keowee measured temperature profiles, 1977 - Station 500 
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Fiqure B-7. Lake Keowee measured temperature profifes, 1977 - Station 500 
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Figure B-8. Lake Keowee measured temperature profiles, 1978 - Station 500 









77 


Figure B 9. Lake Keowee measured temperature profiles, 1971 - Station 501 
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Figure B-10. Lake Keowee measured temperature profiles, 1972 - Station 501 
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Figure B-11. Lake Keowee measured temperature profiles, 1973 - Station 501 
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Figure B- 12. take Keowee measured temperature profiles, 1974 - Station 501 







Figure B- 13. Lake Keowee measured temperature profiles, 1975 - Station 501 
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Figure B- 14, Lake Keowee^ measured temperature profiles, 1976 - Station 501 
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Figure B-1 5. Lake Keowee measured temperature profiles, 1977 — Station 501 
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Figure B-16. Lake Keowee measured temperature profiles, 1978 - Station 501 
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Figure 17. Measured temperature profiles. Lake Keowee, Station 502, 1971 
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Figure B-18.' Lake Keowee measured temperature profiles, 1972 - Station 502 
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Figure B- 19. Lake Keowee measured temperature profiles, 1973 - Station 502 
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Figure B-20. Lake Keowee measured temperature profiles, 1974 - Station 502 
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Figure B-21. Lake Keowee measured temperature profiles, 1975 - Station 502 
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Figure B- 22. Lake Keowee mo»sured temperature profiles, 1976 - Station 502 
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Figure B-23. Lake Keowee measured temperature profiles, 1977 - Station 502 
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Figure B-^24. Lake Keowee measured temperature profiles « 1978 Station 502 
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Figure B-25. Lake Keowee measured temperature profiles, 1971 - Station 503 
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Figure B- 26. Lake Keowee measured temperature profiles, 1972 - Station 503 
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Figure B- 27. Lake Keowee measured temperature profiles, 1973 - Station 503 
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Figure B- 28. Lake Keowee measured temperature profiles, 1974 - Station 503 
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Figure B-29. Lake Keowee measured temperature profiles, 1975 - Station 503 
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Figure B-30. Lake Keowee measured temperature profiles, 1976 - Station 503 
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Figure B-31. Lake Keovree measured temperature profiles, 1977 - Station 503 
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Figure B- 32. Lake Keowee measured temperature profiles, 1978 - Stiition 503 
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Figure B- 33. Lake Koewee measured temperature profiles, 1971 - Station 504 
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Figure B-35. Lake Keowee measured temperature profiles, 1973 - Station 504 
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Figure B-36. Leke Keowee measured temperature profiles, !97<l - Station 504 
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Figure B- 37. Lake Keowee measured temperature profiles, 1975 - Station 509 
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Figure 38. Lake Keowee measured temperature profiles, 1976 - Station 508 
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Figure B-39. Lake Keowee measured temperature profiles, 1977 - Station 504 
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Figure B-40. Lake Keowee measured temperature profiles, 1978 - Station 504 
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Figure B- 41. Lake Keowee measured temperature profiles, 1971 - Station 505 
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Figure B-42. Lake Keowee measured temperature profiles. 1972 - Station 505 
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Figure B-43. Lake Keowee measured temperature profiles « 1973 - Station 505 
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Figure B- 44. Lake Keowee messured temperature profiles, 1974 - Station 505 
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Figure B- *15. Lake Keowee measured temperature profiles, 1975 - Station 505 
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Figure B-46. Lake Keowee measured temperature pirx>files, 1976 - Station 505 
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Figure B-47. Lake Keowee measured temperature profiles, 1977 - Station 505 
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Figure B- 48. Lake Keowee measured temperature profiles, 1978 - Station SOS 
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Figure B-49. Lake Keowee measured temperature profiles, 1971 - Station 506 
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Figure B-50. Lake Keowee measured temperature profiles, 1972 - Station 506 
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Fiqure B-51. Lake Koewee measured temperature profiles, 1973 - Station 506 
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Figure B-52. Lake Keowee measured temperature profiles, 1974 - Station 506 
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Figure 6-54, Lake Koewee measured temperature profiles, 1976 - Station 506 
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Figure B- 55. Lake Keowee measured temperature profiles, 1977 - Station 506 
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Figure B- 56. Lake Keowee measured temperature profiles « 1978 ^ Station 506 




APPENDIX C 


MEASURED TEMPERATURE PROFILES 
(Average of Stations 500 Through 506) 
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Figure C-1. Lake Keowee averaged measured temperature profiles. Stations 501-506, 1971 
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Figure C- 2 . Lake Keowee averaged measured temperature profiles. Stations 501-506, 1972 
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Figure C- 3. Lake Keowee averaged measured temperature profiles. Stations 501-506, 1973 



Fiqure C-4. Lake Keowee measured temperature profiles. Stations 501-506, 1974 





Figure C- 5. Lake Keowee measured temperature profiles. Stations 501-506, 1975 








Figure C-7. Lake Keowee measured temperature profiles. Stations 501'-506, 1977 
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Figure C-8. Lake Keowee measured temperature profiles. Stations 501-506, 1978 




